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Introduction. Abundant evidence is emerging that 
many lavas on Mars were emplaced as slow-moving 
pahoehoe flows (Figure 1). Models for such scenarios 
contrast sharply with those for steep-sloped applica- 
tions where gravity is the dominant force. The mode of 
flow emplacement on low slopes is characterized by 
toe formation and inflation [1-7]. In the latter phase of 
pahoehoe flow emplacement, stagnation, inflation, and 
toe formation are most closely tied to the final topog- 
raphy, dimensions, and morphologic features (Figure 
2). This mode of emplacement is particularly relevant 
to the low slopes of planetary surfaces such as the 
plains of Mars, Io and the Moon. 



Figure 1. Numerous flows on Mars show evidence of 
inflation, similar to terrestrial pahoehoe flows [ 81 . 

The fundamental difficulty in developing a new 
model for pahoehoe lava flows is that random effects 
associated with inflation, internal fluid pressure, and 
crustal strength dominate the emplacement [1-10]. 
Baloga and Glaze [4] developed an initial 2- 
dimensional model based on classical uncorrelated and 
correlated random walks [11]. They showed that for 
simple lobes, the tendency for channel development 



Figure 2. Terrestrial pahoehoe showing randomness in 
toe unit directions and orientations. 


and similar features broadly agreed with field observa- 
tions. Glaze and Baloga [10] presented a more com- 
prehensive model that simulated random transfers of 
individual lava parcels as a function of space and time. 
The output was 3 -dimensional topography showing 
how the lobe thickened and expanded with time sub- 
ject to a variety of factors such as the source geometry 
and flow rate. Of particular interest was rate of infla- 
tion, i.e., the rate of thickening at the expense of ex- 
pansion at the margins. 

T emperature-dependent Simulations. 

The explicit cooling of surface parcels and the in- 
fluence of internal inflation pressure have been added 
to the model of Glaze and Baloga [10]. The three types 
of simulations described here are referred to as “equi- 
probable” (as in [10]), “temperature-dependent”, and 
“pressure-dominated” respectively. All simulations are 
done on a flat pre-existing surface with a constant rate 
of lava supply. Each parcel volume is equal to a typical 
pahoehoe toe (see [9, 10]). 

Ultimately, the mechanical strength of the crust 
controls the movement of lava within the lobe and at 
its margins. The surface temperature is taken here as a 
proxy for the mechanical strength of the crust. The 
surface cooling rates of pahoehoe lavas are well 
known, both from theoretical and empirical studies 
e.g., [2,12-14]. The equiprobable transfer rules in [10] 
have been modified to increase the probability when 
the surface temperature is relatively high and decrease 
transfer probabilities for cooler parcels. 

Mathematically, the probability that basal location 
j at time t is the site of the next lava parcel transfer is 



where T/t) is the temperature of the surface at location 
j and the i summation is taken over all occupied basal 
locations at time t. Once the location for the lava 
transfer has been determined by (1), the algorithm 
treats the four possible directions for the lava transfer 
as equiprobable. The simulations assume that a parcel 
of lava begins to cool when a cell is first occupied and 
thus exposed to the atmosphere. Other parcels trans- 
ferred to that location inflate the volume at that loca- 
tion with the original surface cell continuing to cool 
according to the empirical formula, 

7X0 = -60.81n(0 + 303 (2) 

where t is measured in hours and T is given in °C [2]. 

Internal transfers only inflate the lobe locally and 
leave the pre-existing crust undisturbed and continuing 




to cool. Heat is propagated through the crust very 
slowly [12-14]. Thus, it is assumed that only the sur- 
face parcels (20 cm thick) cool to any significant de- 
gree. Interior parcels remain at a constant temperature 
(~1150°C) until they break out into an unoccupied cell 
at the existing margin of the lobe. The initial surface 
temperature of a breakout at the margin is set to 
1150°C no matter when it occurs in the simulation. 
Subsequently the breakout parcel cools according to 
(2) (Figure 3). 
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Figure 3. Pahoehoe surface cooling curves 
adapted from [2] and [14]. Arrows indicate time step 
nunbers, where the parcel from time step 1 has the 
lowest surface temperature. 

The general topography and morphology of the 
temperature dependent lobes is tightly clustered near 
the origin as in the equiprobable case [10]. Surface 
parcels cool so rapidly that most parcels are on the 
long flat end of the cooling curve and thus have rough- 
ly the same probability of a transfer. 

Pressure-Dependent Simulations. 

When there is an internal lava transfer, the lobe 
thickness at that location increases. This local inflation 
increases the fluid pressure at all connected fluid loca- 
tions within the lobe. The most recent margin breakout 
location would have the weakest confining mechanical 
strength. Consequently, when there are one or more 
internal transfers that cause inflation, the model forces 
the next transfer to occur at the last margin breakout. 

The inclusion of inflation pressure in this way has a 
dramatic effect on the morphology and dimensions of 
the lobe particularly its maximum length (Figure 4). 
Simulated lobes are much more elongated than either 
the equiprobable or temperature-dependent lobes. The 
total volume of all types of simulations has been fixed 
so that the greater length of these lobes results in a 
thinner overall thickness [10]. Despite the changes in 
lobe shape, the lobe volume attributable to inflation is 
essentially the same in all three types of simulations. 

Conclusions: These three types of simulations pro- 

vide a physical basis for interpreting the final stages of 
pahoehoe flow lobes on planetary and terrestrial sur- 


faces. The simulations show that the cooling of indi- 
vidual surface units has little influence on the mor- 
phology and dimensions of a lobe. Specifically, the 
temperature-dependent and equiprobable simulations 
are unlikely to be distinguishable even in the field for 
terrestrial analogs. However, when inflation is present, 
the units take on a marked lobate character that should 
be readily distinguishable in the field or in plan form 
planetary images. It is noteworthy that the rate of infla- 
tion is virtually indistinguishable whether the cooling 
of surface units is ignored, considered, or combined 
with inflation pressure within the lobe. Future studies 
based on HiRISE topography of flows shown in Figure 
1 will attempt to interpret the story of their emplace- 
ment in terms of these simulations. 
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Figure 4. Pressure dependent simulations illustrate 
elongate lobe structure in plan-view. “X” indicates the 
simulation origin. Note warmer surface parcels at the 
flow front. 
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